function y = f(x)
y = x^4-8*x^3-35*x^2+450*x-1001;